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The X-ray Properties of the Energetic Pulsar PSR 
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ABSTRACT 

We present and interpret several new X-ray features of the X-ray pulsar PSR J1838- 
0655. The X-ray data are obtained from the archival data of CHANDRA, RXTE, and 
SUZAKU. We combine all these X-ray data and fit the spectra with different models. 
We find that the joint spectra are difficult to fit with a single power law; a broken 
power-law model with a break at around 6.5 keV can improve the fit significantly. The 
photon index changes from T = 1.0 (below 6.5 keV) to T = 1.5 (above 6.5 keV); this 
indicates a softer spectral behaviour at hard X-rays. The X-ray flux at 2-20 keV is 
found to be 1.6 x 10 -11 ergs cm -2 s _1 . The conversion efficiency from the spin-down 
luminosity is ~ 0.9% at 0.8-10 keV, which is much higher than that (~ 10~ 3 % - 
10~ 4 %) of the pulsars that show similar timing properties. We discuss non-thermal 
radiation mechanisms for the observed high X-ray conversion efficiency and find that 
emission from the magnetosphere of a greatly inclined rotator is the most favorable 
interpretation for the conversion rate and the pulse profiles at X-ray bands. A line 
feature close to 6.65 keV is also detected in the spectra of SUZAKU/X1S; it might be 
the K Q emission of highly ionised Fe surrounding the pulsar. 

Key words: pulsars: general - gamma-rays : observations - X-rays: general - radia- 
tion mechanisms: general - line: identification 



1 INTRODUCTION 

PSR J1838-0655 is an X-ray pulsar that was recently dis- 
covered. This X-ray source was first catalogued by Einstein 
Imaging Proportional Cou nter (IPC) in the Galact ic Plane 
Survey (IE 1835.3-0658; iHertz fc Grindlavl 1 19881 ). ASCA 
found that this source, which was named as AX J1838.0- 
0655 according to the nomenclature of ASCA, is located 
at the sou thern edge of the supernova remnant SNR 
G25.5+0.0 <|Bamba et al.ll2003l ). ISugizaki et al.l (|200ll ) sug- 
gested AX J 1 838. 0- 0655 to be a variab l e sou rce, however, 
iBamba et ail, (|2003f ) and iMalizia et"al] (|2005l ) re-analysed 
the ASCA data and obtained a steady flux for this source. 
AX J1838. 0-0655 was finally identified as an X-ray pulsar 
because of t he discovery of the spin per iod at 70.5 ms using 
RXTE data (|Gotthelf fc Halpern|[20ul ). 

AX J1838. 0-0655 is a bright X-ra y source with a har d 
(r = 0.8 ± 0.4) non-thermal spectrum (IBamba et alj|2003l). 
The observed fluxes of ASCA (0.7-10 keV: IBamba et al.l 
120031 ). Chandra (2-10 keV : iGotthelf fc Halpernl 120081 ) 
and Suzaku (0.4-10 keV; lAnada et all 120091 ) are ~ 
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0.88, 1.1 and 1.32 x 10 11 ergs cm 2 s 1 , respectively. 
The flux in 2-6 keV of the ASCA spectrum (~ 0.4 x 
10 -11 ergs cm~ 2 s _1 ) was different from that of the EX- 
OSAT source GPS 1835-070 (~ 6.3 x 10" 11 ergs cm" 2 s _1 ), 
which was an X-ray source discovered at the positional 
uncertainty of AX J1838.0-0655 with the Medium-Energy 
(ME) proportional counters on EXOSAT (|Warwick et al.l 
Il988l ). However, the variability of the point-source flux be- 
tween the two investigations of EXOSAT and ASCA might 
be caused by the contamination from nearby sources in the 
detection of EXOSAT. The IBIS/ISGRI on board INTE- 
GRA L has also detected this X-ray source a t hard X-rays/7- 
rays. l|Bassani et alj|2004 [Bird et alj|2004 l ■ 

A single absorbed power law usually gave a good fit 
to the X-ray spectra of PSR J1838-0655, but the photon 
indices might vary at different wavebands. A single power 
law with a photon index of T = 1.5 ± 0.2 and a col- 
umn density of N H = (6.7 ± 1.3) x 10 22 cm -2 provided 
a good fit to the composite spec tra from 1 to 300 keV 
generated by ASCA/INTEGRAL (jMalizia et alj|2005l ). At 
low-energy wavebands, the observations of ASCA/G1S and 
Chandra/ ACIS indicated a harder spectrum with F — 0.8 ± 
0.4 (0.7-10 keV) and 0.5 ± 0.2 (2-10 keV) jSugizaki et all 
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l200ll ; iGotthelf fc Halpernl |20b"l ). However, we note that 
the p hoton indices derived from Suzaku/ XIS dAnada et all 
l2009l . T = 1.27 ± 0.11) and Swift/XRT (jLandi et al]|2006l . 
r — 1.86lQ 4g) were significantly softer than the results 
of ASCA/GIS and Chandra/ ACIS at similar wavebands. 
On the other hand, the results of Beppo-SAX/PDS, INTE- 
GRAL/lSGKl and Suzaku/HXD showed much softer spec- 
tra at high-energy wavebands with T arou nd 2.5 to 1.6 
dMalizia et alj|2004 120051; lAnada et al.ll2009h Furt hermore, 
IGotthelf fc Halpernl |200a ) and lAnada et all l|2009h claimed 



that the spectrum of PSR J 1838-0655 became steepened at 
around 8-15 keV; the truth and the origin of the spectral 
break or steepening are still unclear. Besides, we note that 
the detected X-ray emission might contain emission from 
the pulsar wind nebula (P WN) (T = 1.1 - 2.0 in 2-10 keV, 
IGotthelf fc Halpernlliooj ) and/or emission from the super- 
nova remnant. 

The high-energy nature of AX J1838. 0-0 655 has caught 
more attention when lAharonian et al.l (|2005l ) showed a pos- 
sible connection between AX J1838. 0-0 655 and a high- 
energy TeV source, HESS J1837-069. lAharonian et all 
1200611 ) indicated that HESS J1837-069 might be associ- 
ated with supernova remnants or pulsar wind nebulae; and 
this suggested that AX J1838. 0-0655 might be a pulsar can- 
didate before its pulsation was confirmed. Furthermore, the 
non-thermal spectral behaviour with a photon index similar 
to the general pulsars (F ~ 0.5 — 2.5) and the large flux 
ratio between the X-ray an d optical band, -Fb.4-iokev/-FieV 
(|Gotthelf fc Halpernl |200S| ; this value is ~ 100-1000 for the 
Vela pulsar) also gave supporting evidence indicating this 
X-ray source to be a pulsar. 

No radio counterpart of the X-ray pulsar PSR J1838- 
0655 has been detected yet. If HESS J1838-069 is the TeV 
counterpart of the X-ray pulsar, the radio brightness of 
HESS J1838-069 will be extremely faint comparing with 
other known TeV sources. Besides, unlike the Geminga pul- 
sar, PSR J1838-0655 has no detectable 7-rays from hu n- 
dreds MeV to hundreds GeV (|Hartmanll 19991 ; lAbddl2009l ). 

The energy conversion efficiency of PSR J1838-0655 
is very unique comparing with other known X-ray pulsars 
with similar properties. For example, the characteristic age, 
the surface dipole magnetic field, and the spin-down lumi- 
nosity of PSR J1838-0655 are all very similar to those of 
PSRs B0833-45 (the Vela) and B1706-44; however, the en- 
ergy conversion efficiency of PSR J 1838-0655 is about two to 
three order of magnitude higher than those of PSRs B0833- 
45 and B 1706-44 as inferred from the observed X-ray emis- 
sion in 0.8-10 keV. Therefore, PSR J1838-0655 provides 
a unique opportunity to investigate the energy conversion 
mechanism and to understand the nature of the X-ray emis- 
sion of pulsars. 

To understand the puzzling properties of the X-ray 
emission of PSR J1838-0655, we re-examine the X-ray spec- 
tral behaviour of this pulsar using Chandra, RXTE and 
Suzaku archives. We propose three possible mechanisms to 
explain the observed X-ray spectrum and the high energy 
conversion efficiency of PSR J1838-0655. We also marginally 
detect a line feature at around 6.65 keV, which might be Ka 
emission of highly ionised Fe on or around the surface of the 
pulsar. 



2 OBSERVATIONS AND DATA ANALYSIS 

To study the emission mechanisms of PSR J 1838-0658, we 
re-analysed its spectral and temporal behaviours using the 
X-ray archives of Chandra, RXTE and Suzaku. 



2.1 Chandra Data 

The Chandra data were observed on 2006 August 19 us- 
ing the Advanced CCD Imaging Spectrometer operating in 
the Timed Exposure (TE)/VFAINT mode with an expo- 
sure time of 20 ks. The sa me data have been analysed by 
IGotthelf fc Halpernl (|200ST > . The data reduction and spec- 
tral analysis were performed using standard procedures with 
the X-ray packages of CIAO (ver. 3.4 1.1), CALDB 4.1.1 
and X SPEC (ver. 12.4.0). We followed IGotthelf fc Hatoernl 
|2008l ) in selecting the regions for the pulsar, the PWN, and 
the background. The spectra of PSR J1838-0655 and the 
PWN are grouped to have at least 50 counts for each spec- 
tral bin. 



2.2 RXTE Data 

The RXTE data were observed from 2008 February 
17 to March 5 using the Proportional Counter Array 
(PCA; 2-60 keV). The s ame data set has been used by 
IGotthelf fc Halpernl (|2008h to determine the period and pe- 
riod derivative of PSR J1838-0655. The time resolution of 
the data set is very high (1 fis) and can provide a proper 
detection of the pulsed spectra of the pulsars. We performed 
the standard data reduction procedures with the FTOOLS 
(ver 6.5) of HEASoft (ver. 6.5.1) using the latest PCA cal- 
ibration files (20090313). The events that we used for spec- 
tral analysis were restricted to Good Xenon and Good Time 
Interval (GTI). The total exposures of GTIs are 12.9 ks for 
PCU 0, 9.3 ks for PCU 1, 31.5 ks for PCU 2 , 4.1 ks for PCU 
3 and 8.2 ks for PCU 4. We also set the energy boundary of 
photons at 2-20 keV. 

After getting the final products, we performed the so- 
lar system barycentric time correction using "faxbary" at 
(J2000) R.A.=279°.513, decl.=-6°.926 to produce an event 
time list for following analysis. According to th e timing 
ephemeris derived by IGotthelf fc Halpernl l|2008f ). we set 
the period as 70.49824397 ms and the period derivative 
4.925 x 10~ 14 s • s _1 at the epoch of 54522.000696574 MJD 
and folded the RXTE data to produce the pulsed profile of 
PSR J1838-0655 as shown in Fig. [U 

To generate the pulsed spectrum, we assumed the un- 
pulsed emission to be the background. We divided the data 
into "on-peak" and "off-peak" emission depending on their 
pulsed phase. The response of each PCU was produced by 
the Perl script of "pcarsp" and then combined by the task 
of "addrmf" with the weighting depending on the relative 
exposures. The pulsed spectrum was then obtained by sub- 
tracting the off-peak emission from the on-peak emission. 



2.3 Suzaku Data 

The Suzaku observations of PSR J1838-0655 were carried 
out in March of 2007. The archival data contain observa- 
tional results of the X-ray Imaging Spectrometer (XIS) and 



The X-ray Properties of the Energetic Pulsar PSR J 1838-0655 3 




Phase 

Figure 1. (a). The upper panel shows pulse profile in 2-20 kcV 
folded by 32 bins with the period of 70.49824397 ms and the pe- 
riod derivative of 4.925 X 10" 14 using 8 RXTE/PCA data from 
2008 Feb. 17 to Mar. 5. The time zero (to) of the profile was 
set at the epoch of 54522.000696574 MJD. (b).The lower panel 
shows the pulse profile in 10-60 kcV folded by 16 bins with the 
period of 70.496725 ms using Suza/cu/HXD-PIN data of 2007 
Mar. 5. The time zero of the profile was set at the start of GTI 
(54164.53394270 MJD). Both (a), and (b). are divided into two 
groups to represent the pulsed phase and unpulscd phase. We de- 
termined the pulsed phase (on-peak emission) as 0.625 cycle. The 
pulsed phase of (a), starts from 1st to 20th bin of 32 bins and 
it of (b). starts from 4th to 14th bin of 16 bins. The plot of the 
counts are normalized to the average photons in each bin. 



the non-imaging Hard X-ray Detector (HX D). The data sets 
have been analysed bv lAnada et~aT] (|2009l ). 



2.3.1 XIS 

The Suzaku/XIS observations of PSR J1838-0655 have a 
total exposure of 42.2 ks. The data were observed in the 
normal mode without window option and the pixels on the 
CCD were read out every 8 sec. The 42.2 ks exposure of XIS 
observation was divided into two editing modes of 5x5 and 
3x3 with 14.9 ks and 27.3 ks respectively. 

The selection of source area for the Suz aku/XIS data 
in our analysis is slightly different from that of lAnada et al.l 
l|2009D . Because PSR J1838-0655 is contaminated by some 
nearby sources, which can margi nally be resolved in the XI S 
image (e.g. GPSR5 25.252-0. 139; iGotthelf fc Halpernl200sl ). 
we extracted the spectra of PSR J1838-0655 only within l'.5 
circular regions centered at (J2000) R.A -18 h 38 m 03M3, 



decl.=— 06°55'33".4 instead of the regular adoption of 3' cir- 
cular region for a point source. Our selection still contains ~ 
70% of the total energy for a point source. We extracted the 
background of the pulsar from a concentric annulus of radii 
2'. 5 < r < 3'. 5 centered at the pulsar position and set the 
energy boundary to be 0.5-10 keV. The extracted spectra 
comprise 2935 counts for XIS 0, 2180 counts for XIS 1 and 
3427 counts for XIS 3 after subtracting the background. 

The XIS data reduction and spectral analysis were ex- 
ecuted using XSELECT (ver. 2.4a) of HEASoft (ver. 6.5.1) 
and XSPEC (ver. 12.4.0) with the latest Suzaku/XIS cal- 
ibration files (20090203). We generated the response ma- 
trix (rmf) and auxiliary response (arf) files with the HEA- 
Soft command of "xisrmfgen" and "xissimarfgen" . The spec- 
tra are produced with each channel containing at least 50 
counts. A cross-calibration term was included to correct the 
difference among three XIS spectra at 0.5-10 keV. 



2.3.2 HXD 

The Suzaku/HXD observations of PSR J1838-0655 have a 
total exposure of 37.7 ks and offer a time resolution of 
6 /is for pulsed detection. However, the spectrum of the 
Suzaku/GSO (Gadolinium Silicate phoswich counters) data 
are not good enough for statistical analysis, we thus only re- 
duced and analysed the pulsed spectral data obtained from 
the HXD-PIN (Positive Intrinsic Negative silicon diodes) 
data. We performed the data reduction and spectral analy- 
sis procedures of the PIN data using XSELECT (ver. 2.4a), 
FTOOLS (ver 6.5) of HEASoft (ver. 6.5.1) and XSPEC (ver. 
12.4.0) with the latest HXD calibration Files (20090203). 
The events were restricted to the effective energy range (10— 
60 keV) of the PIN detector. We also applied solar system 
barycentric time correction with the task "aebarycen" at 
(J2000) R.A -279°. 513, decl.=-6°.926 in order to produce 
and analyse the pulsation of PSR J 1838-0655. Based on the 
timing ephemeris derived bv lGotthelf fc Halpernl (|2008h . we 
set a trial period of 70.496725 ms at epoch MJD 54164.98049 
(the mid-point of this Suzaku/HXD observation) . 

The hard X-ray spectrum obtained by subtracting the 
non-X-ray background and cosmic X-ray background was 
divided into two groups of phase cycles to represent the "on- 
peak" and "off-peak" emission. The pulsed spectral data 
were obtained by subtracting the off-peak emission from the 
on-peak emission and rebinned to ensure that the photons in 
each spectral channel are larger than 50 after subtracting the 
un-pulsed background. The pulsed photons of PSR J1838- 
0655 above 50 keV are very few in the data and are thus 
ignored in spectral analysis. 



3 RESULTS 

Table [1] shows the results of our spectral analysis for Chan- 
dra, RXTE, and Suzaku observations. We also show the re- 
sults of previous analysis in Table 1 for comparisons. Our 
results are consistent with previous analyses; however, some 
significant discoveries are found in our work. First, we iden- 
tify a spectral break at 6.5T1.0 keV, which is roughly consis - 
tent with the previous claim of IGotthelf fc Halpernl (|2008T) 
and lAnada et alj (|2009l) that the spectrum becomes steepen 
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Table 1. Single Power-law spectra for AX J1838.0-0655/PSR J1838-0655 detected by different X-ray observations. 



Instrument 


Energy range 


Photon index 


Column density 


Unabsorbed flux 


References 




(keV) 


(T) 


(10 22 cm" 2 ) 


(10 — 11 ergs cm -2 s —1 ) 




Einstein/IPC 


0.2-3.5 








0.014 a 


Hertz & Grindlav (19881 


ASCA/GIS 


0.7-10.0 


0.8 (0.4-1.2) 


4.0 (2.8-5.7) 


1.1 


Sugizaki et al. (20011 


Beppo-SAX/PBS 


20-100 


2.5 (2.2-2.8) 




3.6 


Malizia et al. (20041 


INTEGRAL/lSGRl 


20-300 


1.66 (1.43-1.89) 




9.0 


Malizia et al. (20051 


Swift/XKT 


0.2-8.5 


1.86 (1.4-2.55) 


5.54 (3.73-8.37) 


1.13 6 


Landi et al. (20061 


Chandra/ ACIS 


2-10 


0.5 (0.3-0.7) 


4.5 (3.7-5.2) 


0.88 


Gotthelf & Halpern (20081 






0.5 (0.2-0.7) 


4.1 (2.7-5.5) 


0.86 


This work 


RXTE/PCA C 


2-20 


1.2 (1.1-1.3) 


4.5 (fixed) 


0.9 b 


Gotthelf & Halpern (20081 






1.3 (1.2-1.4) 


4.1 (fixed) 


0.92 b 


This work 


Suzaku/XIS 


0.4-10 


1.27 (1.16-1.38) 


5.4 (4.9-5.9) 


1.32 


Anada et al. (20091 




0.5-10 


1.37 (1.28-1.50) 


6.0 (5.6-6.5) 


1.26 


This work 


Suzaku/HXD C 


12-50 


2.0 (1.1-3.0) 




1.84 


Anada et al. (20091 




10-50 


2.0 (0.9-3.3) 


6.7 (fixed) 


1.9 


This work 



Notes: The errors are set in 90% confidence interval throughout 

a) X-ray flux in IPC counts s . b) Indicated in 2—10 keV. c) Only the pulsed spectra were calculated. 



at around 8-15 keV. Second, we find a line feature around 
6.4 keV in the spectra of Suzaku/yilS. 

We have applied a simultaneous fit to a compos- 
ite spectrum of PSR J1838-0655. We included the spec- 
tra of ACIS / Chandra, PC A/ RXTE and HXD / Suzaku; all 
these spectra are believed to be dominated by the pulsed 
phases. To account for the cross-calibration mismatch be- 
tween each instrument, we also introduced a constant to 
this fit. In order to compare with the the joint spectrum 
of ASC 'A/INTEGRAL data (|Malizia et al.ll2005h , we have 
fixed the column density to 6.7 x 10 22 cm -2 in our analy- 
sis. A single absorbed power-law model for the composite 
spectrum have a power-law index of Y = 1.2 ± 0.1 with 
xt = 1-35 for 91 dof. The statistics does not improve even 
we set the absorbed column density as the free parameter 
to fit the spectra. 

A broken power-law model can obviously improve the 
fit of the joint spectrum and the additional power-law com- 
ponent is significant at more than 99% via a F-test. The 
power-law indices change from I.O^q'j to 1.5l '2 with the 
broken energy at 6.51*'q KeV. The statistics have signifi- 
cantly improved to = 1-20 for 89 dof. We note that a 
simi lar spectral break h ave been discovered in the Crab pul- 
sar (|Kuiper et al.ll200ll 1. 

To check whether the broken power-law is caused by 
the instrumental responses of different observations, we have 
also applied the power-law models to different energy bands 
of the RXTE data to examine the variation of the pho- 
ton indices of the pulsed spectrum. Originally, no signifi- 
cant variation of the photon index was found between the 
energy domains less than 10 keV and larger than 10 keV 
(r ~ 1.0 — 1.1). However, we note that both photon indices 
are flatter than that of the total spectrum. We thus examine 
the photon index of the spectrum at the energy of 5-15 keV 
and find the photon index T ~ 1.4. This suggests that there 
is a spectral break at around 10 keV, and the break is not 
caused by the uncertainty of the response functions among 
the different instruments. 

We also note that the RXTE spectrum is very dif- 
ferent from the Chandra one, which shows a much flatter 
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Figure 2. Cross-calibration fit of the joint spectra (Blue: Chan- 
dra/ ACIS; Red: RXTE/PCA and Green: Suzaku/HXD-PIN ; The 
vertical axis of the RXTE spectrum is referred to only one PCU). 
We applied a broken power-law with a change of the slope from 
1.0 to 1.5 and a break energy at 6.5 keV to get the best spectral 
fit. The bottom panel shows the residuals in terms of as with 
error bars of size one. 

photon index (F ~ 0.5). We note that the time resolution 
for the Chandra TE mode data is only about 3.2 sec, and 
it is not adequate to get the pulsed spectrum from tim- 
ing analysis. However, if the un-pulsed emission is mainly 
contributed by the surrounding ionised cloud, the pulsed 
spectrum can roughly be obtained by considering the spec- 
trum of its PWN as the background of the pulsar. The 
PWN spectrum can be well fitted with a power law with 
a photon index of F ~ 1.5(1.1 — 2.0), and the flux of the 
PWN is ~ 1.6 x 10~ 12 ergs cm -2 s -1 . The flux ratio between 
PSR J1838-0655 and its PWN is > 5.3 in our calculation, 
and this ratio is at th e high end of the distributio n of PWNe 
observed by Chandra (jKargaltsev fc Pavlov||2008l l . This sug- 
gests that the observed pointed-source spectrum of Chandra 
is mainly from the pulsar even we do not have an appropriate 
time selection. This suggests that the un-appropriate time 
selection can not be the cause of the different the photon 
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Figure 3. Spectra of 3 XIS detectors from 2-10 keV (Blue: XIS 0; 
Green: XIS 1 and Red: XIS 3). Each spectral channel has at least 
30 counts. In the simultaneous fitting for the spectra of two FI XIS 
detectors, we proposed the absorbed power-law and Lorentz line 
profile to fit the data sets and the line emission was marginally 
detected close to 6.4 keV. The bottom panel shows the residuals 
in terms of as with error bars of size one. 

indexes between RXTE and Chandra spectra. The origin 
of the difference is still unclear; but since the Chandra and 
RXTE observations were carried out at different epochs, the 
difference might indicate that the photon index could vary 
with time if the difference is not caused by the uncertainty 
of the response functions between the different instruments. 

A very interesting new discovery from the Suzaku obser- 
vations is a marginal line feature close to 6.4 keV detected 
in the XIS spectra of 3x3 mode (Fig. [3}. This line feature is 
only detectable in the FI detectors (XIS and XIS 3) and 
can not be seen in the spectrum of the BI detector (XIS 1); 
this might be caused by the fact that the effective area of the 
BI chip above 4 keV is smaller than that of the FI chips and 
thus the XIS 1 is much less efficient in detecting line emis- 
sion at higher energy. Detailed analysis show that this line 
feature is at around 6.65 keV (EW = I50tl1° eV ) and have 
a significance of ~ 4cr. This line feature might be the K a 
emission of highly ionised Fe. We do not find similar emis- 
sion feature in the spectra of the surrounding background. 
This excludes the possibility that this line feature may be 
caused by the background contamination, such as the intense 
diffuse Galactic ridge emission. This line feature might be 
originated from the surrounding ionised gas of the pulsar. 



4 DISCUSSION 

We have presented the results of temporal and spectral anal- 
yses for the X-ray emission from PSR J 1838-0655 observed 
by CHANDRA, RXTE and SUZAKU. According to Fig. [T] 
the structure of the pulse profile of Suzaku/PIN observation 
is very similar to that of RXTE/PCA. There is some dis - 
similarity of the pulse profile derived bv lAnada et al.l (|2009l ) 
and ours. This may be caused by the different selection of 
the zero epoch and the phase bins. We also divided the pho- 
tons of RXTE observations into 2-10 keV and 10-20 keV to 
plot the folded light curves (shown in Fig. [4| , and we did 
not find obvious difference from the profile structures. 




Phase 

Figure 4. Pulse profile of RXTE/PCA observations folded at 
different energy bands. All the input parameters are correspond- 
ing to Fig. [T] The upper panel shows the folded light curve of 
PSR J1838-0655 at 2-10 keV with ~ 70.5 ms and the lower panel 
shows similar folded profile for the same pulsar at 10—20 keV with 
the same period. The plot of the counts are also normalized to 
the average photons in each bin. 

Our spectral analysis indicates that there is a spectral 
break in the X-ray spectra of the pulsar. This might be 
simi lar to the spectral break discovered in the Crab pul- 
sar l|Kuiper et alj|200lf ). With non-thermal emission model, 
the best fit photon indices are Fi ~ 1 and F2 ~ 1.5 
with a break energy E c ~ 6.5 keV and with the flux of 
/ ~ 1.6 X 10~ n ergs cm" 2 s" 1 in 2-20 keV. The results pro- 
vide more complete information of the spectral behav i our in 
additio n to the studies do ne by iGotthelf fe Halpernl (|200Sf ) 
and bv lAnada et al.l (|2009l ). At a distance of d ~ 6.6 kpc, the 
efficiency of the energy conversion from the spin-down lumi- 
nosity, L p — 5.5 x 10 36 ergs/s, to the X-ray emissions is about 
0.9%, where we have used -Lo.s-iokev ~ 5.0 x 10 34 erg/s es- 
timated from our data analysis. 

If the observed X-ray emission is caused by a non- 
thermal process, the emission might come from the mag- 
netosphere around a particle acceleration region, which is 
called a "gap" . Comparing the efficiency of the energy con- 
version from the spin-down energy to the non-thermal X- 
ray emission with other known non-th ermal X-ray emitting 
pulsar s listed in the tables 1. and 2. of iKargaltsev fe Pavlov! 
(|2008f ). we find that the effi ciency of PSR J1838-0 6 55, w hich 
was listed in the table 3 of IKargaltsev fe Pavlovl (|2008l ) for 
pulsars without known periods, is the third largest and 
is larger than that of the Crab pulsar. One might expect 
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Figure 5. Optical depths from the pair-creation of the 7-ray 
photons as a function of energy. The solid line and the dashed 
line represent the optical depths for PSR J1838-0655 and the 
Crab pulsar, respectively. 



that the energy-conversion efficiency of PSR J 1838-0655 
should be similar to that of the 7-ray pulsars, PSRs B0833- 
44 (the Vela) and B1706-44, because these three pulsars 
have similar pulsar properties (e.g. the rotation periods and 
the surface magnetic fields). However, the efficiencies of the 
Vela pulsar and PSR B1706-44 are about 5 x 1CT 4 % and 
3 x 1CP 3 %, respectively, and are much smaller than that of 
PSR J1838-0655. In fact, the efficiency of PSR J1838-0655 is 
in the range of younger 7-ray pulsars, such as the Crab pul- 
sar and PSR B1509-58, which have Lo.8-iokev/£p ~ 0.34%, 
and ~ 0.61%, respectively. Therefore, PSR J1838-0655 pro- 
vides a unique opportunity to investigate the non-thermal 
emitting process in the pulsar magnetosphere and the con- 
nection between the X-ray and 7-ray emission. 

According to t he outer magnetospheric emission model 
(|Cheng et alJl200Ch . the high efficiency of the Crab pulsar 
can be explained as a result of the X-ray photons emitted 
via the synchrotron radiation of the secondary electrons and 
positrons produced in the pair-creation process of the pri- 
mary 7-rays, which are produced in the gap. In the gap of 
the Crab pulsar, about 1% of the spin-down energy is con- 
verted to the 7-ray emission. Considering the optical depth 
of the pair-creation of the 7-rays using the X-ray photon 
number density outside the gap inferred from the Crab ob- 
servations, we find that most of the primary 7-ray photons 
above 1 GeV are absorbed by the background X-ray photons 
before escaping from the magnetosphere (Fig. [5}. Because 
there are abundant new born electron-positron pairs, which 
will emit X-ray photons via the synchrotron radiation, the 
X-ray emission will dominate the observed energy spectral 
distribution of the Crab pulsar and the resultant conversion 
efficiency in the X-ray bands can be as large as a few point 
percent. 

The above emission model predicts the photon indices 
to be Pi ~ 1.5 and T2 ~ 2 with a break at E c ~ 10 keV 



l| Takata et al. I [20071 ). Roughly speaking, the predicted spec- 
tral property is also similar to the observed properties of the 
X-ray emission of PSR J1838-0655, which shows a spectral 
break around 6.5 keV and has a high energy conversion effi- 
ciency. However, we find that the X-ray emission model for 
the Crab pulsar can not be applied to the case of PSR J 1838- 
0655 because the optical depth of the pair-creation from the 
7-rays is much smaller than unity. Using the observed X-ray 
properties in this study, the energy distribution of the X-ray 
photon number density in the magnetosphere of PSR J 1838- 
0655 is 

( {E x /6.5keY)- 1 , 

dN o fi m22j2 -3 -1 I E x < 6.5keV . . 

— =2.8x10 d 6 . 6kpc cm erg I ( ^ /g ^-1* (1) 

{E x > 6.5keV 

We have used = —-rS-i^Y 2 , where is the 

dEx c dEx v 2 ' ' dEx 

observed number flux obtained in this paper, and Ri c is the 
light cylinder. Fig. [5]compares the optical depths of the pair- 
creation for the 7-ray photons of the PSR J1838-0655 (solid 
line) and of the Crab pulsar (dashed line) before escaping 
the magnetosphere. We see that the optical depth is much 
smaller than unity for the GeV 7-ray photons emitted in the 
gap of PSR J1838-0655, while it is about unity for those of 
the Crab pulsar. So the absorption of the 7-ray photons in 
the magnetosphere of PSR J1838-0655 is very weak, and 
only a few secondary pairs are created by the pair-creation 
process with the background X-ray fields. As a result, the 
synchrotron emission of the secondary pairs will not be able 
to explain the flux of the observed X-ray emission. On this 
ground, we conclude that the X-ray emission mechanism of 
the Crab pulsar can not apply to PSR J1838-0655, although 
the energy-conversion efficiencies of both sources are similar. 

A greatly inclined rotator, with the inclination angle 
a ~ 90°, is another possibility to explain the high energy 
conversion in the X-ray emission of PSR J 1838-0655 and 
the different efficiencies between PSR J1838-0655 and the 
Vela pulsar. With the outer magnetospheric model, the X- 
ray emission from the Vela pulsar can be explained by the 
synchrotron emission of the secondary pairs. The secondary 
pairs are produced by the pair-creation process between the 
7-rays and the background X-ray field and/or by the mag- 
netic pair-creation of the inwardly propagating 7-rays, which 
pass through the region near the stellar sur face. The outer 
magne tospheric emission model studied by IZhang fc Jiand 
|2006t ) argued that the efficiency of the X-ray emission in- 
creases with the inclination angle. This is caused by the 
fact that (i) the strong acceleration of the outer gap model 
extends between the null charge surface of the Goldreich- 
Julian charge density and the light cylinder, and (ii) the 
null charge surface shifts toward the stellar surface with in- 
creasing inclination angles from 0° to 90°. For the greater 
inclined rotator, therefore, more 7-rays emitted inward di- 
rections pass through vicinity of the stellar surface and are 
absorbed by the strong magnetic field around the stellar sur- 
face. The created secondary pairs can emit X-ray photons 
via the synchrotron radiation. 

The magnetic pair-creation condition for a photon with 
the energy E~, may be written as ( — ■ \ ) -§- sin 6t, ~ 0.2 

(Muslimov fc Hardind I2003F ). where B c = 4.413 x 10 13 
Gauss and 9^ is the angle between the magnetic field lines 
and the propagating direction of 7-rays. Most of the pairs 
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Figure 6. Cascade of the magnetic pair-creation process. The 
vertical axis represents the energy of the emitted photons, and 
the abscissa axis represents the generation number. The dashed, 
solid and dotted lines represent the pair-cascade initiated by the 
photons with 10 GeV, 1 GeV and 0.1 GeV, respectively, emitted 
inwardly at the radial distance 7*0 = 0.055 in units of the light 
radius. The photon energies of the generation number i represent 
the typical energy of the synchrotron photons, and the values of 
Ti in the figure represent the radial distances to the pair-creation 
positions of i-th generation. 



in the gap are cr eated around the null charge surface 
ijTakata et al.l[20o3 ). The distance from the stellar surface 
to the null charge surface of the last-open field line is typ- 



ically r, 

(3 tan a + ^/9tan 2 a + S) 



sin 2 (6 



and tan ( 



a)sin#( c , where tan(9 n = 

_ (3+^9+8 tan 2 q) 



2 ^v, c — -tana 

with 9 n being the angle of the null charge surface measured 
from the north pole and 9i c being the angle to the point at 
which the last-open field line is tangent to the light cylinder. 
If we use the radial distance r n as the emission point, the 
closest distance from the stellar surface to the trajectory of 
the inwardly propagating 7-rays emitted on the null charge 
surface is r m ;„ ~ r n sin#j,, where 6 P is the angle of the null 
charge surface measured from rotation axis. If we estimate 
the minimum inclination angle, above which most of 1 GeV 
7-rays are absorbed by the strong magnetic field, with the 
magnetic pair-creation condition, we find a > 60°. In real- 
ity, the pairs are created above the last-open field lines; for 
example, if we set the emission point with 2r„, we would 
obtain a > 70°. 

The absorbed 7-ray photons are converted into electron- 
positron pairs. The new born pairs lose their energy rapidly 
via the synchrotron radiation; and the synchrotron photons 
would be converted into electron-positron pairs via the mag- 
netic pair-creation process, if their energies are large enough 
for the magnetic pair-creation process. The typical energy 
of the synchrotron photons emitted by the first pairs are 
Ei ~ 0.075-E 7 , where E-y is the energy of 7-ray photons. 
The typical energy of the synchrotron photons emitted by 
the zth generated pair can be written as Ei ~ (0.075)*B 7 . 



If the magnetic fields on the path of the photons are 
strong enough, the pair-creation cascade will continue un- 
til the typical energy of the synchrotron photons becomes 
about 1 MeV. Let us take the inclination angle of 80° and 
a dipole magnetic field as an example. The radial distance 
to the null charge surface on the last-open field line in the 
magnetic meridional plane is r n ~ 0.028i?; c (Figure [6]). We 
compute the pair-creation cascade initiated by 1 GeV pho- 
tons emitted inwardly along the magnetic field line at the 
radial distance r ~ 0.055i?; c and on the null charge sur- 
face in the magnetic meridional plane. As summarized in 
Figure [S] the first pair-creation occurs at r = 0.014i?i c and 
the typical energy of synchrotron photons emitted by the 
first generated pairs can be written as E\ ~ 75 MeV. The 
sequence of the pair-creation cascade stops after producing 
the 3rd generation pairs at r/ ~ 0.004J?; C . The 3rd genera- 
tion pairs are produced with a Lorentz factor of 7 ~ 6 and 
emit synchrotron photons of ~ 500 keV. The spectrum of 
the emission from the pairs extends with a photon index of 
1.5 - 2 (Cheng & Zhang 1996) from ~ 500 keV to a break 
energy E m ~ 14 keV, which is the characteristic energy of 
the synchrotron photons from the electrons/positrons with a 
Lorentz factor 7 = 1/ sin 6b ~ 1. Below the break energy, the 
spectrum will have a photon index of F ~ 2/3, which repre- 
sents the spectral slope of the synchrotron radiation below 
its characteristic energy. We find this model can explain the 
observed properties in the X-ray spectra of PSR J1838-0655. 

Therefore, if the X-rays from PSR J1838-0655 are emit- 
ted in the magnetosphere, an inclination angle close to a ~ 
90° would be preferred in order to explain the energy con- 
version rate at X-ray bands. For the Vela pulsar, the energy 
conversion rate can also be e xplained by a = 30° with the 
magnetic pair-cascade model (|Zhang fc JiandboOfjf ) . On this 
ground, it is possible to explain the difference of the energy 
conversion efficiencies of the Vela pulsar and PSR B1838- 
0655 by the difference of the inclination angles. The greatly 
inclined rotator model may also predict the absence of the 
radio detection. The copious pairs produced by the mag- 
netic pair-creation near the pulsar can easily supply the 
Goldreich- Julian number density. Therefore, if the magnetic 
pair-creation cascade initiated by the 7-ray from the outer 
gap develops in the polar cap accelerator, the electric field 
parallel to the magnetic field in the polar cap accelerator will 
be screened by the new born pairs, indicating no coherent 
radio emission in the magnetosphere. 

We note that most of 7-ray photons emitted outward 
can escape from the magnetosphere without pair creations. 
The outwardly migrating particles are accelerated using a 
full potential drop along the field line before escaping the 
gap from the outer boundary located around the light cylin- 
der; on the other hand, the inwardly migrating particles are 
accelerated with a small part of the potential drop between 
the inner boundary and the pair-creation position around 
the null charge surface. This implies the luminosity of the 
outward 7-rays are one or two-order larger than that of the 
inward 7-rays. If the observed X-rays are originated from 
the inward 7-rays, the outward 7-rays would be observed 
with a flux similar to or greater than the X-ray flux, which 
is fx ~ 10 -11 ergs cm -2 s _1 . However, we note that the ob- 
served 7-ray flux also depends on the viewing angle of the 
observers. Therefore, future Fermi observations could thus 
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provide more information necessary to investigate the X-ray 
emission process with a greatly inclined rotator model. 



5 CONCLUSIONS 

We have re-analysed the X-ray data of CHANDRA, RXTE, 
and SUZAKU to investigate the possible X-ray emission 
mechanism for the X-ray pulsar PSR J1838-0655. We find 
that a broken power-law with the photon index varying from 
1.0 to 1.5 and with a break energy of ~ 6.5 keV can provide 
a good fit for the composite X-ray spectrum. According to 
our results, we also found no significant energy dependence 
of the pulse profiles. 

We have discovered a marginal Fe K a line feature in 
the Suzaku/XIS data of PSR J1838-0655. A similar line 
feature was also detecte d in the radio pulsa r PSR J 1420- 
6048 /AX J1420. 1-6049 ((Roberts et alJl200ll ) in ASCA ob- 
servations. Both pulsar are found to be associated with high- 
energ y TeV sources (i e HES S J1420-607 for PSR J1420- 
6048; lAharonian et al.1 l2006al ): however, PSR J1420-6048 
not only has a radio counterpart but also a hundred-MeV 
7-ray one (3EG J1420-6048). These results might indicate 
that the Fe K a emission is associated with the TeV sources; 
however, the relation between the line emission and the TeV 
sources is still not clear. 

We have proposed and investigated the possible mech- 
anism for the observed X-ray emission of PSR J1838-0655. 
The model for a greatly inclined rotator is the most plausible 
scenario to explain the high efficiency of the energy conver- 
sion from the spin-down luminosity to X-ray emission. How- 
ever, it is expected to accompany with some level of 7-ray 
radiation. Because no unidentified EGRET source exists in 
the region around PSR J1838-0655, it requires a deeper ob- 
servation in the 7-ray bands (i.g., Fermi telescope) to solve 
the mystery of the X-ray emission from this pulsar. This will 
also provide important clues to understand how the X-rays 
and the 7-rays in the magnetosphere of PSR J1838-0655 are 
connected. 



ACKNOWLEDGMENTS 

The authors thank an anonymous referee for his/her help- 
ful comments. We also thank Dr. Albert Kong for the 
careful review of the manuscript and thank Drs. Hsiang- 
Kuang Chang, Kwong-Sang Cheng and Ronald Taam for 
fruitful discussion. This research has made use of the data 
obtained through the High Energy Astrophysics Science 
Archive Research Center Online Service, provided by the 
NASA/Goddard Space Flight Center. This work was par- 
tially supported by the National Science Council through 
grants NSC 98-2811-M-008-044. CYH acknowledges sup- 
port from the National Science Council through grants 
NSC 96-2112-M-008-017-MY3 and NSC 95-2923-M-008-001- 
MY3. JT was supported by the Theoretical Institute for Ad- 
vanced Reserch in Astrophysics (TIARA), operated under 
Academia Sinica and the National Science Council Excel- 
lence Projects program in Taiwan through grant NSC 96- 
2752-M-007-007-PAE. 



REFERENCES 

Abdo A. A., 2009, ArXiv e-prints 

Aharonian F., Akhperjanian A. G., Aye K.-M., Bazer- 

Bachi A. R., et al. 2005, Science, 307, 1938 
Aharonian F., Akhperjanian A. G., Bazer-Bachi A. R., 

Beilicke M., et al. 2006a, A&A, 456, 245 
Aharonian F., Akhperjanian A. C, Bazer-Bachi A. R., 

Beilicke M., et al. 2006b, ApJ, 636, 777 
Anada T., Ebisawa K., Dotani T., Bamba A., 2009, PASJ, 

61, 183 

Bamba A., Ueno M., Koyama K., Yamauchi S., 2003, ApJ, 
589, 253 

Bassani L., Malizia A., Stephen J. B., Bazzano A., Ubertini 
P., Barlow E. J., Bird A. J., Dean A. J., Hill A. B., Shaw 
S. E., Walter R., Bodaghee A., Belanger G., Lebrun F., 
Terrier R., Much R., 2004, The Astronomer's Telegram, 
232, 1 

Bird A. J., Barlow E. J., Bassani L., Bazzano A., Bodaghee 
A., Capitanio F., Cocchi M., Del Santo M., Dean A. J., 
Hill A. B., Lebrun F., Malaguti C, Malizia A., Much R., 
Shaw S. E., Stephen J. B., Terrier R., Ubertini P., Walter 
R., 2004, ApJ, 607, L33 

Cheng K. S., Ruderman M., Zhang L., 2000, ApJ, 537, 964 

Gotthelf E. V., Halpern J. P., 2008, ApJ, 681, 515 

Hartman R. C. e., 1999, ApJS, 123, 79 

Hertz P., Grindlay J. E., 1988, AJ, 96, 233 

Kargaltsev O., Pavlov G. G., 2008, in Bassa C, Wang Z., 
Cumming A., Kaspi V. M., eds, 40 Years of Pulsars: Mil- 
lisecond Pulsars, Magnetars and More Vol. 983 of Amer- 
ican Institute of Physics Conference Series, Pulsar Wind 
Nebulae in the Chandra Era. pp 171-185 

Kuiper L., Hermsen W., Cusumano G., Diehl R., 
Schonfelder V., Strong A., Bennett K., McConnell M. L., 
2001, A&A, 378, 918 

Landi R., Bassani L., Malizia A., Masetti N., Stephen J. B., 
Bazzano A., Ubertini P., Bird A. J., Dean A. J., 2006, 
ApJ, 651, 190 

Malizia A., Bassani L., Landi R., Molina M., Stephen J. B., 
Gianotti F., Schiavone F., Barlow E. J., Bazzano A., Bird 
A. J., Capitanio F., Dean A. J., Del Santo M., Lebrun F., 
Renaud M., Shaw S. E., Terrier R., Ubertini P., Walter 
R., 2004, in Schoenfelder V., Lichti G., Winkler C, eds, 
5th INTEGRAL Workshop on the INTEGRAL Universe 
Vol. 552 of ESA Special Publication, From BeppoSAX to 
INTEGRAL: PDS Observations of Hard X-Ray Sources 
Detected in the IBIS Survey, p. 161 

Malizia A., Bassani L., Stephen J. B., Bazzano A., Ubertini 
P., Bird A. J., Dean A. J., Sguera V., Renaud M., Walter 
R., Gianotti F., 2005, ApJ, 630, L157 

Muslimov A. G., Harding A. K., 2003, ApJ, 588, 430 

Roberts M. S. E., Romani R. W., Johnston S., 2001, ApJL, 
561, L187 

Sugizaki M., Mitsuda K., Kaneda H., Matsuzaki K., Ya- 
mauchi S., Koyama K., 2001, ApJS, 134, 77 
Takata J., Chang H.-K., Cheng K. S., 2007, ApJ, 656, 1044 
Takata J., Shibata S., Hirotani K., 2004, MNRAS, 354, 
1120 

Warwick R. S., Norton A. J., Turner M. J. L., Watson 

M. G., Willingale R., 1988, MNRAS, 232, 551 
Zhang L., Jiang Z. J., 2006, A&A, 454, 537 



